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Abstract 

A collision singularity in the N-body problem has the property 
that all position vectors have finite limits and at least one mutual 
distance approaches zero as time approaches the instant of the 
singularity. General properties of a coxiision singularity are inves- 
tigated and it is fo^nd that they are analogous to those of a collision 
of all bodies at the same point and time. A more detailed study of 
a simple type of such a singularity, namely that of several simulta- 
neous binary collisions, is presented. 
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Introduction 

Little is known about the singularities of the equations of motion 
of the N-body problem of celestial mechanics, particularly if complex 
values of the independent and dependent variables are admitted* Only 
the three simplest cases have received considerable attention, as 
manifested by the bibliography: the binary collision in the three -body 
problem (e* g. , Sundman [4]) and in the N-body problem (e. g. , 

Wintner [2]), the triple collision in the three-body problem (e. g* , 
Sundman [1]» Siegel [!])» and the collision of all N bodies at the 
same point and instant (e* g. , Block [2], Wintner [2]); all of the above 
studies are restricted to real values of all variables. For a study of 
the binary collision in the N-body problem for complex values of the 
variables, see Sperling [3]. 

It is presently unknown which types of real singularities (in the 
sense of complex analysis as well as with regard to their dynamical 
interpretation) are possible in the N-body problem; the only general 
result is a remarkable statement by von Zeipel [2] (cf. 205), 
asserting that a singularity is a collision singularity (cf* 203) if the 
mutual distances remain bounded as time approaches the instant of 
the singularity. **Collision singularities** are such that the N bodies 
separate into distinct ‘‘clusters** with all bodies in each cluster 
colliding at the instant of the singularity at a well-defined point, 
while the clusters remain apart from each other. Since von Zeipel’s 
statement seems to be virtually unknown and particularly in view of 
Wintner*s [2, p. 431] remark: ** A note of H. von Zeipel indicates a 
consideration to the effect that, if U becomes infinite when t tends 
to a finite value, then J must tend to infinity, unless all bodies tend 
to definite limiting positions. But it seems to be hard to fill in the 
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gaps " , we v/ill subsequently present a detailed proof, following 
essentially von Zeipel's ideas, Chazy's [6] attempt of a proof is 
only a brief sketch and must be considered at least vague. The 
introduction of " cluster coordinates " (cf, 103-104) enables us to 
prove in a simple manner that the existence of finite limits of the 
mutual distances implies the existence of the limits of the positions 
and relative positions as t approaches a finite value (cf. 216); 

Wintner [2, p. 327] mentions this as an undecided problem (cf, von 
Zeipel [2, p. 4]). 

Ic will be shown that the known important properties of a 
collision of all bodies at the same instant hold, with obvious modifi- 
cations, also for a collision singularity, the main result probably 
being that the bodies of each cluster are positioned close to a central 
configuration when the time is close to the instant of the singularity. 

A detailed study of the case of several simultaneous binary 
collisions concludes the paper. This case has already been considered 
by Lahaye [3]; his treatment, however, is at least incomplete: his 
assumption, that the asymptotic behavior of a perturbed two -body prob- 
lem equals that of the two -body problem as the instant of a binary 
collision is approached, is not justified without proof, and he does 
not show that his solutions are the only ones. The main result of the 
investigation is that this singularity gives rise to an algebraic branch 
point of the second order on the real time axis, the same as for one 
binary collision; in particular, the motion can be analytically continued 
in the real beyond this singularity. 
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While writing this paper, I was informed (an abstract in the 
Notices Amer. Math, Soc, 15 (1968) 6^1 and private communication) 
that H. Pollard has worked in the same field. Details of his work 
have not been r ailable to me at the time of the completion of this 
paper. 


1 

1 
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On the notation 


Let t^^ ^±00 (later t^^ - 0) be the instant of the (collision) singularity; 
all limits are taken as t -* t^' ( t < t'^^ ). 


C , 

c 

denote 

positive 

constants 

1 

t^ , bounded on an interval 

b , 

b 

denote 

functions 

of 

t or s = 






[to. t 

with suitable < t^s 

e , 


denote 

functions 

of 

t or s. 

with e 0 , e 0 as t t* 


The above symbols have fixed meaning, at least within the given 
context, if they are indexed; without indexes they are unspecified, 
i. e. , they may denote different constants or functions from one 

occurrence to the next . For example, in 

t 

^b(T)dT = b(t) and in f(t) = b(t) + tb(t) the b*s 
0 

occurring in the same equation may denote different 
functions. 

1 

Ti.e absolute value of a vector z is z ( = Izl = +(zz)^), 

^ denotes asymptotic behavior as t t^M= 0) or s — 0 : f ^ g is 

equivalent to f = g + o(g) . 

Differentiation with respect to time t is denoted by ' : df/dt = T . 

Differentiation with respect to s is denoted by ' : df/ds f' . 

The gradient of a scalar f with respect to a vector ^ is denoted 
by grad^f = 

Subscript following semicolon: 0: some fixed value; 1, 2, 3: the three 
components of a vector. 
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100 Equations and iiitegrais of motion 

For reference we list here the equations and integrals of motion 
and some related concepts and equations, in the two coordinate systems 
used later. The geometric basis is the 3 -dimensional euclidean space. 

101 (Cartesian) inertial coordinates, the center of mass resting at the 
origin. N particles (^bodies) with masses m^ . > 0 are moving 
according to Newton's law of gravitation. Let y be the gravitational 
constant, the position vector from the origin to the body , 

and 

z., = zv - z* (101. 1) 

the relative position vector from m^ to mj^ . 

Set ]vj 


M = 


.m. 


( 101 . 2 ) 


k=l 


102 Equations of motion 


nv " Sj 

N 

-vE 

-jk 

Zk--v2^m. 

m. 

J 3 

j = l l^k-Sj' 

j=l 

2jk 

ji^k 

j^k 


5U 




k = 1 -- N 


( 102 . 1 ) 


potential energy 


^ ^ 2-^ I y., - z ■ I 


-vEE-4 


m. m, 


l^j<k^N 


jk 


( 102 . 2 ) 


kinetic energy 
N 

V = m- 

j = l 


N 


j=l 


(102. 3a) 
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N N 

V = “ 777 " 7 / m.m., z'-i,z\ , 

4M / ^ /— < j k’-JK-JK. 

j=l k=l 

polar inertia momentum 

N N 

J = 7 m.z-z. =7 m-z? 

^ •j'-j-j Lu J j 

j=i j=i 


'em 

j=l k=l 


£ 2 ] "‘j”k5jk2jk =1^X11 ”"j 


m, z^ 


l<j<ksN 


k jk ’ 


angular momentum 
N 


K 


E m. z- X z'. = 

J"J •“j 


N N 

1 \ V 

2M ”^j”^k5jk^Sjk 

j=l k=l 


conservation of energy (energy theorem) 
V + U = h ( = constant of energy) 
conservation of angular momentum 


K - constant , 

center of mass rests at the origin 
N 

L m- z. = 0 , 

j = l 

Lagrange -Jacobi equation 

j'' = -2 U + 4h = 2 V + 2h 


the positions Zt are expressed i’l terms of the relative 
N“^ 

by V' 


j = l 


Jk 


(102. 3b) 

(102. 4a) 
(102. 4b) 

(102.5) 

( 102 . 6 ) 

(102.7) 

( 102 . 8 ) 

(102.9) 

positions 


( 102 . 10 ) 
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103 "Cluster coordinates"* Partition the set of all N bodies into L 

nonempty disjunct subsets, called "clusters" in the following. Let , 

1 < N < N , be the number of bodies in the jJ-th cluster. Obviously, 

+ ... + Nl =N. (Observe that subscripts are separated by commas 

in cluster coordinates and that = k-I-- ). Notation: 

JL 


m = mass of the k-th body in the jJ-th cluster , 
z. ^ ~ position vector from the origin to m. , , 

= position vector from the origin to the center of mass of the 
£-th cluster, 

r. , - position vector from the center of mass of the jJ-th cluster 

to m , . 
k 



center of mass 
of jK-th cluster 


Obviously , 


z 

«4* 


I, k 




(103. 1) 


define M , and Z . by 
N. 


N. 


At A t 

E m and M.Z = / m ,z , ; (103. 

i, k £mmi £,k-ji,k 


k=l 


2 ) 


k=l 
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then N 


t 


m c . 1 =0 
i— ^ SL, ki.^. k 

k=l 


ar 


N 


M 


1 

i£i,k ^iJ j^S^ik 

j=l 


(103. 3) 


(103.4) 


104 Equations of motion 


•''ZZ"' 




z , - z, 
k H 


k 


X> K j 13 

X = l„=l l£i,k-Sx,»l 

&, k^\, H 


(104. 1) 


L ”x ^ 


M Z'' = -Y 


z. , - z. 

"•A.k <“X, K 


ZZZ^i.k-^x.K, 

X=lH = lk=l '«5x,k’5 J 


(104.2) 


«^,k ' «X,H ” ijt.k "■£x,h 


c''' 

ijJ.k 


■vZ 


ni , 


£j^,k £x,j 


j 


J ?^k 


>■ ”x 

zz 

X -1 K-l 
X 


N. 


m, 


X, H 


(104. 3a) 


' >"x £je,k ~.kx,K _.j_ Y“ 

z. - z, +c. ,. -c. I" ‘"‘iz, -z, +c 


Z - Z + C - f 

***X -AiU ix,v 

m. ■"■■ " ■ “ 

4^ 


*Z ’^X wi,k »^Xik' 


. r 
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or 




1 u . 


^ — + 




. ^ 


(104. 3b) 


Define the potential energy of the £-th cluster by 


N N 
I i 


m . .m 


U 


V'V' 

& - -^2^Lj~ — : — : 

k < j 




(104.4) 


its kinetic energy by 
N 

a 

“j.kfLkC.k 

k- . 

its (total) energy by 


(104.5) 


. . U 

I 0. i 


(observe that ^ is not 

I 


(104.6) 


constant in general) , 


its polar inertia momentum by 
N. 


a 

SL JL^ I, k:^j£, kX£, 
k=l 


(104.7a/ 


N N 

i a 

* (104.7b) 

^ j=l k=l 
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and its angular momentum by 
<«£ ^ ^ ^£j0,k 


(104.8) 


define further the kinetic energy ^f the centers of mass of all 
clusters by 


V' - i/ M z'z' 

L I 


(104.9) 


their polar inertia momentum by 
L 




MfZ.Z, 


(104. 10) 


and their angular momentum by 




Z„ X z, 

I ■'•I 


(104. 11) 


write 




+ u+ = U* + 


(104. 12) 


observing that U'*’ consists exactly of those terms of U which contain 
the difference of coordinates of masses in different clusters. 


Using (103 3;, it is easy to show that 


V = ^ + V° = V* + v° 


(1C4. 13) 


+ J° = J* + J° 


(104. 14) 
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and 


L 

K = ^ 
Jj=l 



+ 


K* + K° 


L 

h = + V° + U-^ = h* + V’ + U* 

£ = 1 


( 104 . 15 ) 


( 104 . 16 ) 


200 Classification and properties of real singularities 

As mentioned before, ^ve restrict ourselves (necessarily) to real 
values of the independent variable time t and the dependent variables, 
the positions and velocities; consider onl/ real analytic continuation of 
the solution of the equations of motion and exclude the two cases that 
t = (i)®, or that the finite t is an accumulation point of singularities 
from both sides: then, if there is a singiilarity at t^' , this point is 
the right or left endpoint of an open interval on which the 

solution is holomorphic. Without restriction of generality, we will assume 
that t^' is the right endpoint of such an intervals 
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201 For definiteness, we refer to the equations of motion (102, 1), but 
analogous statements obviously hold for the equations in the form 
(104.2) , (104.3). Let t^ , >2.k*o* ^k*o finite and real such that 

I z, - ^0 for jj^k; the existence theorem for analytic 

' ‘*k;o >o ^ 

differentia.1 equations there exists an interval 1^ “ ^ol ^ ^ 
a unique holomorphic s<-lu.tion Zj^(t) ; a k = 1 --N, of (1C2. 1) 

satisfying noteworthy (and soon to 

be used) that an estimate of 6 can be made involving only 

minlz, - Z:„l and max| z\ I , but not z, (cf. Siegel "3T, 

'**k;o •“JjO' 'tMk;o' *“k;o ' & ^ j» 

pp. 23 -25). Continue thij solution holomorphically from t^ for 

increasing t ; the finite point t* is a singularity, if the solution is 

holomorphic on [t^»t*[ , but not on [t^»t*]. 


202 Theorem a): t* is a singularity iff liminf(minz., ) = 0 as t -» t* , 


b); t* is a singularity iff lim(minz.,) = 0 as t-»t*. 


Proof: 


1, If t* is a singularity then lim inf(min z^^) = 0: assume contrariwise 
that lim inf(min z^^^) s dg > 0 . 

Then there exist numbers tj and d^ , tg s tj < t* and 0< d^ £ d^ , 
such that ^jk ^ infer that | | = ^ ^ dj"® 

and, using the equations of motion (102.1), ISk^ ^ • From 


V 

5k 5k^^ ) + (t - tj ) z^(tj ) + f (t - t) z^(t) dT 


it follows that z, and z' have finite limits z,'" and z'**' with 

«k *«k “k **k 

z* = I zf - zf| s d, for ji^k as t t* . The existence theorem 
jk -*k *“j' ~ 
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(cf. 201), applied to initial values very close to t* , z!'' , states 

that the solution remains holomorphic at t* , contrary to the assumption. 


2. If lim inf(min = 0 , then t^' is a singularity: if at least one 

is unbounded as t t^' , then there is a singularity at t^' . Assume 
now that all remain bounded; then, by integration, all are 


bounded. In 


dz, . dz . 

j = l k=l k;j k;j 


are bounded by (102. 1) and the 


assumption; thus, and therefore U are bounded. This is incon- 
sistent with the fact that liminf(minz ) = 0 implies limsupU = — oo , 

3. If is a singularity, then lim(minz.- ) = 0: assume contrariwise 

that lim sup(min z^j^) ^ d^ > 0. There exists a sequence t^ t^ 
such that z^j^(t^) ^d^>0, — U^dg, and by the energy theorem 

V < d 3 , implying |zMt )| ^ d^ . Apply the existence theorem (cf. 201) 

k V 

to t , z, (t ) , z^ (t ) as initial values; if t is sufficiently close 
V V '-k' V V 

to t^' , t* is contained in the interval 1 1 - t | < 6 where the solution 
remains holomorphic, yielding a contradiction (observe that 6 is a 
function of d^ and d^ only, but not of t^). 


4. If lim(minz ) = 0 , then t^ is a singularity: this is implied by 2. 


Observe that a) does not require the energy theorem (102.6) for its 
proof. 
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203 Definition: A singularity t* is called a collision singularity if all 
position vectors z, have finite limits as t t^' . 

Hence, by 202 theorem, at least one of the mutual distances z., 
converges to zero as t -• t^' for this type of singularity. 

According to the original definition, the partition of the bodies into 

clusters (cf. 103) is quite arbitrary; in the following this partition is 

always used in conjunction with a collision singularity and assumed to 

be consistent with it, i. e, : two bodies m. and m^^ are in the same 

cluster iff z._ is very small as t t^' = instant of the collision 
jk 

sing'ilarity, or in terms of cluster coordinates: 


all 







are very small compared with 
on the considered time interval ♦ 


(203. 1) 


204 Lemma: Let the z^ = holomorphic on tQ < t < t^'. 

Then lim J ^ <» exists as t t^' , 

Proof (cf, IVintner [2] p, 327): 

Let t^' be a singularity, otherwise the statement is trivial; then 
lim(minz., ) = 0 , hence lim(— U) = +<» , From = 4h - 2U it follows 

that j'"' is ultimately positive; hence, is ultimately increasing and 
therefore ultimately does not change sign, and we conclude that J is 
ultimately monotonic (increasing or decreasing), implying the existence 
of the limit of J . Since J ^ 0 , limJ is <» or a nonnegative 


number 
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205 Theorem (von Zeipel): If there is a singularity at t^ and 

lim J < ® as t t* , then t* is a collision singularity. 

Proof: Assume that the solution of the N-body problem is holomorphic 
on "Iq, t-[ , to < t^' . 

206 Let first be limJ = 0 : then Zj^ 0 and 0 by (102,4) , 

implying that 0 and C ; all positions and all relative 

positions z_ have the finite limit zero and we have a collision 
v*jk 

singularity. All bodies collide at time t* at the center of mass. 


207 Now let 0 < lim J = J(t^') < ; we partition the set of all mutual 

distances into four disjoint subsets as follows: 

first type: limz^j^> 0 exists; 

second type: limz,^ = 0 exists; 

jk 

third type: limz.. does not exist and liminfz,, > 0 ; 

jk jk 

fourth type: limz.. does not exist and lim inf z., = 0 , 

jk )k 

lim supz.. > 0 . 
jk 

We show in the following that 0 < lim J < ® implies the existence 

and finiteness of all limz.. . In a first case (cf. 208), we assume 

jk 

that there are no mutual distances of the fourth type and prove our 
assertion directly; in a second case (cf. 209-215), we assume that 
there is at least one distance of the fourth type and lead this 
assumption to a contradiction. 216 lemma then concludes the proof. 


208 First case: there are no mutual distances of the fourth type. 


There exist a t^Q , t^ ^ t^^ < t* , and two constants 
0 < a^ < ag , such that on [tQo,t5l«[ 


and ag , 


for all mutual distances of the first and third type: ^jk ^ ' 

and for all mutual distances of the second type: z^^ < a^ . 
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We consider now the solution on [tQQ,t’^[ • Introduce cluster coordinates 

such that two bodies are in the same cluster iff for their mutual 

distance z.. < a, holds. 

Jk ^ 

Note: Because the mutual distances z., are continuous on 

Jk 

and either z <3.^ ^ u ^ ^ » ^ke 

Separation of the bodies into clusters is unique on this 
interval ; i. e. , each body stays in one and the same 
cluster. 

All mutual distances of bodies in the same cluster, say the £-th, are 
by the definition of the cluster of the second type, i. e. , If .l-^O 

as t -► t^' « By (103. 4) this implies C , 0 , hence C u ® > i. e. , 

Jo » M , K •*♦ 

al? limC exist (and are zero). 

vi«J 09 K 

Consider the equation of motion (104.2) for the center of mass 

of the £*th cluster; on the right of this equation only relative positions 
and mutual distances of bodies in different clusters occur, i. e. , the 
mutual distances are of the first or thix’d type; hence |z^ , ~ z I >Up , 

l(z^ , - z, )|z . - z, and Iz'' I < yklag"^. This 

“A, k -X,H “‘£,k ^ '**A T -s 

together with 




*'0 0 


implies that limZ, exists and is finite as t t^ . We conclude that 

limz^ - =limZ^ *f limC, , $ limz., , limz- , and limz., all exist 
k u,Z, k -^jk ’ k * jk 

and are finite; i. e. , in this first case all mutual distances have 

finite lim?^s. 
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209 Second case: There is at least one mutual distance of the fourth 
type. The existence of lim J = J(t*) > 0 implies that for every e > 0 
there exists a t^_,_ , tg 5 t^^, < u-*- , such that 

|j(t) - J(t*)| < e on [t^^,t>:‘[ . (209. 1) 


Assume without restriction that e <'|-J(t*) ; then J__j_>§J(t*) 
on [t^^,it*[ , and with (102.4) it follows that 



and by (102. 10) 

k = 


^ rto^, t*[ 

3-4 t-^O** * 

--N, on [to^, t#[ . 


(209.2) 
(209. 3) 

(209.4) 


210 Define n(r» d) as the number of mutual distances z., for 

which z., (t) s d holds; consider n(T, d) on the set 
jh 

{ T, d : t^^., ST<t*, 0<dsa3 }; -|N(N-1) being the total number of 


mutual distances, n(T, d) satisfies 


1 5 n(T» d) s ^N(N-l) on ti/is ' .u : 
the first inequality holding because of (209.2). 
Define 

3l(t, d) = g.l.b.6(t, d) 

S(t, d) = { n(r, d) : T€[t, t*[ , t^^, s t < t* , 


( 210 . 1 ) 


( 210 . 2 ) 

0 < d 5 ag , 
t fixed , d fixed ] ; 
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since each n{r» d) is a natural number satisfying (210.1), their 
greatest lower bound exists and is the minimum of all competing n(T, d). 
5(J{t, d) satisfies 

1 5 2Ut, d) <: |N(N-1) - 1 on { t, d : to., ^ t < t* , 0 < d ^ a^ ) , 

(210.3) 

the second inequality holding because of lim(min 2 ^j^) = 0 . 31 (t, d) is 

thus a number such that at every point of the interval [t, t*[ at least 
51(t, d) mutual distances are 2 : d . Furthermore, 

3l(ta, dg) 2 3Ut, d) for t 5 tg < t* , 0 < dg s d ; (210.4) 

for, 6(tg,dg) is a subset of 6(t, d) and hence g. 1. b.6(tg, dg)> 
g, 1. b.6(t, d). 

Define 

31* = 1. u. b. { 3Ut, d) ; s t < t* , 0 < d ^ ag } ; (210. 5a) 

since each 31 (t, d) is a natural number satisfying (210.3) , their least 
upper bound exists and is the maximum of all competing 31(t, d) . 

Let t„ and d„ , t^^ ^ t„ < t* , 0 < d^ s ag , be such that 

9l*=3Ut.,d„) ; (210.5b) 

by (210.4) we have 

31* = 3l(t„ , d„ ) = 3l(t, d) for t. ^ t < t* , 0 < d < d„ . (210.6) 

Note: There are always at least 3l(t, d) mutxial distances with 
z > d on [t, t*[ , but the individual mutual distances satisfying 
this may change from one instant to another* 
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211 We define now the constants used in the subsequent analysis and 
further restrict in this context the interval on which we consider the 
solution. First, choose t^ » t^^ ^ t^. . , < , and a^ > 0 such that 

for all mutual distances of the first and of the third type (i. e, , those 

with lim inf z.. > 0 ) 

Jk 

Zjk>as on » (211,1) 

and that for all mutual distances of the fourth type 


lim sup > a^ on • (211.2) 

By (210.6) djj can be replaced by a smaller positive number with- 
out changing ; replace d^ by d^,^ (subsequently defined). Then 


at every point of least 

31^' mutual distances are ^ d^, , 


Set 


= (io-° i 

' min ' 




= (:o 


-Q-10 i 




) 


® d.. 


and determin ? the natural number Q such that 

d^ < min(d„ , ag ) and < 10 ; 

set further 




. a. 

\ ^ JL min ^ 



1 

s 


9 


(211.3) 


(211.4) 


(211.5) 


( 211 . 6 ) 


(211. 7) 


a, = ag + - t*| Ynd"*" 


( 211 . 8 ) 
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ar- 


2 M^a.; 


+ 


V M 



e - 



10 


-2Q-10 




( 211 . 9 ) 

( 211 . 10 ) 


Now r. t.,, , t^,, < t* , can be de* rmined such that 


(1) t.,, s max(tQ^,^, , t„ , t* - eas ) , 


(2) all inutual distances of the second type satisfy 

z.k < . and 

(3) Ij - J(t*)l < e with e from (211.10) . 


( 211 . 11 ) 


We show: there exist t^ , tg , 1^ , t^,^ < tj^j^ < tg < < t* , such that 

(a) every mutual distance of the fourth type is > ag at least 

once on [ti.,. tg] ; 

(b) exactly ^N(N-l) - Dl* mutual distances are < e^., 

on [tg.t^] . 

(a); Select a t^.,^ on and consider the mutual distances 

z of the fourth type for t increasing from tj^,^ ; since they 

p^ "* 

satisfy limsupz > cig (cf. (211.2)) , for each of them there is an 
P9 

instant t on [t, ,,t*[ such that z (t ) > a» . Let U. , 

pq IV pq pq “ 


tj^. ^ tg^ < t* , be such that on each of them there is 

such an instant t where z (t ) > au ( L. = max(t ) will do) . 

pq pq pq pq 






(b): Consider ail mutual distances on • Assume that at 

every point of this interval more than mutual distances 

are 2 * i. e. that n(r, e^) > 31* , t€ [ta 5 j.*t*[ ; then 3l{t^^, ejj,)>3f*, 

a contradiction to the maximum property of 31* (cf. (210.6), (211.6), 

(211.11)). .fence n(T^, e*) S 31* for infinitely many in C*a^»^*C» 

and since 31(1^^^, d^) = 31* by the choice of t^^. , e^j., and d^, we 

have also n(T, e^) s n(r. d^j.) s 31* for t € [ta:j£* follows that 

n(r^» Cj*.) = 31* , i. e. exactly ^N(N-l) -31* mutual distances are <e^ 

at ; take one of these for , then there is a nearby ^ >tg 

such that by their continuity the number of mutual distances being 

<e^ on [^,^1 still equals ^N(N-1) - SS* . 
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212 We have shown that ^N(N-l) -31^' mutual distances are < 

on rtgjtg] , while all other 31^' mutual distances are > d,, on this 
interval. Set 




€ 


❖ 




on [tg.tg] } 


* 


( 212 . 1 ) 

®d, " ^ ■ ^jk- [ta.ta]} ; 

because the mutual distances ate continuous functions of time, there 
is no exchange of elements between 6 and S • We know that all 

mutual distances of the first and of thci third type are > d^^ on [t 2 ,t*[ 


and therefore in S. (cf. (211. 1)) , a?ad all of the second type < 

d.v 

■V 

on and therefore in S (cf. (211.11)) . We prove that at 

n* 

least one mutual distance of the fourth type must be in 6 , i. e. , 

n" 

< on [tgfts] • assume contrariwise that all (and at least one) 
mutual distances of the fourth type are in 6 . None of the mutual 

distances in S (all of them are of the second type) can increase 

beyond on [t^jt^f , i. e. , no mutual distance can change from 

6 to S , and therefore none of the 31^' mutual distances in 6 , 
d, d. 


can decrease below d. , because then less than 31^ mutual distances 
would be a d» : while bv the definition of 31^^ always at least 51^' 
mutual distances are 2 d,, on . This last result is contrary 

to the foregoing assumption that there exists a mutual distance of the 
fourth type in , because this must necessarily decrease below d^, 

on rtq.t^'r • Thus we ha,ve shown that there is a mutual distance z 

L'a » L pq 

of the fourth type in © • 
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It follows that there is a tg, , I 3 < < t* , such that z = d . 

Consider the mutual distances in S for t increasing from tg ; 
it follows now that 

there is a t^ , 1 ^ < t^ < t=p , such that 

il) at least one mutual distance of S equals d at t. 

and * ( 212 . 2 ) 

(21 all mutual distances of 6 are < d^, on [t 2 ,t 4 ^[ • 

e* 

Considering the mutual distances in S for t decreasing from tg , 

6 <*> 

we find analogously (cf, 211 (a)) that 

there is a t^ , < t^ < t^ , such that 

(1) at least one mutxaal distance of S equals d^,, at t^ 

and (212.3) 

( 2 ) all mutual distances of 0 are < J • 


Observe that 

all mutual distances of remain 2 d^j, on [tj.t^] . (212.4) 

V 

Introduce cluster coordinates such that two bodies are in the same 
cluster iff their mutual distance is < , i. e* , it is in 0 (, and 

Sou 

in different clusters iff their mutual distance is ^ , i. e. , it is 

in ) on [tg.l^ ] . 

By (104.7) and (10414) we have 

J* < on (212.5) 

and 3 

m . 

df at t. and at t^ ; (212.6) 


it follows that has a minimum on » 

2 


m 


- J* . > 
max mm 


min 

2M 


( 1 - 10 “^“^°) 


and (cf. (211.5)) 
on [t^ , t^ ] . (212, 7) 
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213 J° has a maximum on : with J = J - J* we find 

J° > J(t«) - e - el = J(t«) (1 - 2- on 

2 " (213.1) 


and 


m 


J < J(t^) + e - 


df = J(t*) (1 - 10'^ + 10’^^'^°) 

4 M 


at tj and at t^ ; (213.2) 


the right side of (213.1) is greater than the right side of (213.2) 
(recall that Q :> 1 , natural). 


214 Let a maximum of J , existing by 213 , be at tj , t^ < t^ < t^ 

(214. i) 


J°'(t,) = 0 


) S 0 


We estimate on [t^.t^] : by (104,10) 

L L 

= 2y^ UZ\Z\ + 2y^ M Z Z'' ; 

x=i i=l 


(214.2) 


all z, are boiuided (cf. (209.4)), therefore, all Z are bounded 

k -V 

by (103.2) : 


<^4 


on [tj.t^] 


Estimate of Z^ (cf. (212.4)) : 

I (z , - z )lz , - z. |'*^lsd,,^ for 
k ~X,K^ k **X, k' '' 

and from the equation of motion (104.2) for Z^ we find 

lz"l<vMd;' 

Estimate of z' : 

at tc , (214,1) and (214.2) imply that 


(214.3) 


(214. 4) 


(214. 5) 
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M 


X = 1 

hence (cf. (211.7)) 


X. , L 

z' M z' z\ < ry^M Z Z'' 


< M Y M d^. ' 


Q(t=)l <(v 


z'.(t.)l <(yM^M/ . = a, 

' ' ^ mm * 


')'■ 


6 » 


and from 


= S;^) + / 


it follows that (cf, (211.8)) 




-s 


+ 1 1 - t* I yMd"® = a, on [t^ , t^ ] 


thus finally (cf. (211,9)) 


|j°"|< 2 M^a 7 +YMa 4 _djj,®^ = ag on [t^.t^] 


215 Write (cf. (209. 1)) 

j = j(t*) + , I I < e 


and 


J° = J°(tc ) + £ 


2 * 


I 

/ 


.ovv 


£2 =/ (t - t)J (r)dT 
tc 


^2 I I ^ “ ^4 I ^8 ^ I I ^8 ^ ^ (cf. (211. ] 


Hence , 


J« = j . j° = J(t«) . J°(tg ) + £j - fig 


(214.6) 


; (214.7) 
(214.8) 
(215. 1) 


(215.2) 

1)) . 


9 


(215.3) 
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and 

- J'-^ . < 4e = 2 M € f 

max mm 


2. lo-Q-lO df 


(215.4) 


this being inconsistent with the inequality (212.7). 

Thus, the assumption that there is a mutual distance of the fourth 
type is false, and taking 208 into account, we deduce that all mutual 
distances have finite limits. The proof is concluded by the following 


Lemma: Let the N-body motion be holomorphic on [t^.t^'C , 

to < t^ > and let all mutual distances have finite 

limits as t t^' • Then all relative positions 

and all positions have finite limits as t t* • 


I. e. , at the motion remains holomorphic, or there occurs a 

collision singularity. The proof is analogous to that in 208, observing 
that the set of mutual distances of the third type is empty. 
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300 Properties of a collision singularity 

301 This part generalizes results which have been established for a 
collision of all bodies ("simultaneous collision") , to each cluster of 
the configuration developing prior to a collision singularity. The 
principal result is that, for t close to t* , the configuration of 
each cluster is close to a central configuration. Our proofs are 
similar to those of Wintner, with some extensions and modifications. 


302 Observe that a collision singularity at t* is characterized by: 
at least one mutual distance converges to zero, and J > 0 remains 
bounded and has a limit as t -* t* . Without restriction of generality 
we may assume that t* = 0 and that the approach is from t < 0 , 
Also, without mentioning it always, we assume that t < 0 is very 
close to the instant of collision singularity t = 0 , such that for the 
cluster coordinates holds; 

^Z,k ^ very small on [t, 0[ , and 
I I > a certain positive constant on [t, 0[ , JL^X , 

It follows also that Z^' is bounded (cf. (104.2)), consequently lim Z^ 
and lim Z^ , and finally lim^^'‘ exist as t -* 0 » as seen 
(104.2). 


303 Differential equation and ultimate behavior of J* : 

Using (102. 9): 4h - 2U , (104. 12): U ^ U* + U+ , and (104. 14): 

J = J* + J° , we find 

J*'' = -2U* + 4h -2U+ -1°"' ; (303. 1) 


both U* and 1°'’’* (cf. 214) are bounded on the considered interval, 
and we write 




(303. 2) 
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From J*" = -2U* + b and -U* -* <» as t-»0 we conclude that: 

J*'' +0O , 

J*' is monotonically increasing and ultimately does not (303. 3) 
change sign, and 

J* is ultimately monotonically decreasing. 

The last statement follows from the facts that J* is ultimately 
monotonic (since does not change sign) and J* > C , J* -♦ 0. 

304 Estimate of K* : 

The differential equation (104.3) for Co v and the definition (104.8) 


K yield 








(304. 1) 

k=l 

k=l 



Write, expanding in a binomial series, 

'Si - + il.r -ix. 1 + '^txrsl • 

where ^i\xs absolutely convergent series, each term of which 

contains a or a component of a Cp . as a factor; it is e:isily 

seen that the terms in P, ^ not containing such a factor cancel out, 
i. e. , P. ^ = 0 if all Cb , = 0 , and we conclude that each term of 
P^ ^ contains a factor of the typec mentioned above. 

From the definition (104.7) of J* and (104kl4) 

Cp^, <cJ*^ , Cp,,Cr,s<cJ* , (304.3) 

and the same inequalities for the components of Cp , » hence, 

•H ^ 

1 k' 1 < cJ* and k' = J* b . (304. 4) 

and observing that (cf. (303. 3)) 

T € [t,0[ . 


Integrating from t to t = 0 
J* = J*(t) = maxJ»i‘(T) , 


(304. 5) 
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we find 


and 


K, = K,.^ + tJ*:=b 
‘^1,0 


K*!* - K*!* + t J*‘* b , 

‘-VO 


= K-:*® + tJ-b 
o 


305 ^ d) : 


Using on 


L N, 


£=1 k=l 
the inequality 

(Z^v\)= 


with 


h.k ■ *=v =”i.kC'i,k • 

we find 
and similarly 

Using (ab)® + (axb)® = a®b® , addition yields 

• M ^ tU» 

ij«v3 + J*TTm^ , .c] , = 2J*V* ,, 

4 ^4- A, k.«i, kiii, k (305. 1) 

(cf. (104. 5) , (104. 13)) . 


(304. 6) 
(304.7) 


Substitute for K*® from (304.7) and from (104.13): V = V* + V° , 
from (104. 14); J = J* + J“ , and (102. 9): J*'' = 2h + 2V ; then 

i j»NS + K*s s J*( -2h + J®'' - 2V° - tb) , 

4 o ' ' ’ 
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and since -2h + - 2V*^ - to is bounded on [t, 0[ , we finally get 


5 J* ( J*" + d ) 

4 o 


(305.2) 


306 The inequality (305. 2) implies that K* = ^ and that 
lira J=I- ^ J > 0 exists : 


define 


i 

Q = dJ*® + 


(— + K*^)/j 

4 o 


&2 


(306. 1) 


then 

Q'= |J*"^ + J*'' - ( Y + K*2)/J* 

and since 0 , we see iiiat and Q' are of the same sign, 

implying that Q ultimately caimot increase, since ultimately 

decreases; hence, Q has a limit Q as t-»0, 0^Q<«>, and 

o o 

tiiis implies because of -♦ 0 that 

-i 

=0 and (0^) limJ* ^ (<®) exists • (306.2) 

o 



307 ^ c> 0 : 


= -2U^^ + b implies 

= -2J*^U* + J*2 b ; 

(307. 1) 

now 



, -U* > c(min^^^ , 

therefore > c> 0 . The second 

term in (307. 1) goes to 

zero, and the desired inequality 


1 

J*2 ^ c > 0 

(307.2) 


follows 
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i 1 -i. 

308 c > 0 and — 0=:='^^ ^ J* ( J*'' + d) imply limj* 0 : 

deiine (cf. (306. 1)) 

Q = d + -^ J=r2j*> a ; 

4 

then 

1 1 

QJ*2 = dJ* + — J*'*^ 

t 

= dj* + I f J*' J-'Mt , since J* - 0 by (306. 2) . 

J 

0 

Substitute from (307. 2), then 

t 

1 P r «i 1 

QJ*2 s dj* + ^ J* = dJ’f' + cj*^ 

0 

and this implies 

Q s dJ*^ + c 

i. e., 

1 -— 

LxiiQ = — limJ* 0 . (308. 1) 

4 


303 From (308, 1) asymptotic expressions for J* and as t -• 0 

are easily derived: 

(308. 1) is equivalent to 

= 4ii- + e , e = e(t) - 0 ; (309. 1) 


separating and integrating we find 



Substitution into (309. 1) yields 


(309.2) 
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J*' = + t® e 

or (309. 3) 


310 ultimately: 

differentiation of (303. 1) yields 

J*‘M = _2U*' - ZU"^' - J"'" 


(310. 1) 


by (104,10) J®''' is a linear combination of terms z''‘ and Z Z"''', 

Z , Z' , and Z'^ are bounded, and by the equation of motion (104.2) 

for Z , Z'''‘ is a linear combination of the components of the f'' ^ 

^ Z ^ Z Pi 

with bounded coefficients, plus a bounded function. Since I c'‘ I < cV*® 

U*Pj ^ 

and this inequality holds also for the components of , we get 

utP» 9 


I < c 


(310. 2) 


Similarly, U"*" is a linear combination of the components of the ^ 
with bounded coefficients, plus a bounded function, hence 


U*'* I < c 


'310. 3) 


is a linear combination of terms (r -r )(r^ -c'' )\c -C 

^ P» Q ^ $ s q *-*p» s Mtpf q I 

this term's absolute value bein^ |c' - c' 1 I C_ -C I”® ; since 

“ bip, q «p, s' ’ip,q Sp, s' 


1-3 


1 C - C „ 1 < c U*® , we have 

*p, q »♦?, s 

|U*'| < cV*^U«2 


(310. 4) 


By means of the energy equation U* + V* = h - U* - V® and the 
Lagrange-Jacobi equation (303. 2): = -2U* + b , both V* and U*!* 
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can be replaced by J*", and we finally get the desired inequality 


I J*"' I < c 


(310.5) 


311 Asymptotic expression for J*" . We prove J*® : 

5 j*(j*" + d\ (cf, (305.2), (306.2)) implies 
4 

i 1 

liminfJ*^J*" > — limJ* ; (311.1) 

it remains to be shov/n that 
1 

lira sup . (311.2) 

Introduce 

F = ; 

then (311.3) 


and as consequence of (309. 3) we have 
F-^.12u=^^2t 

1 

Consider now F^ (^ 0 obviously) instead of 

i 

liminfF^= 31iminf(J^' ^ )(J=J^^ ) 

1 

= 12^=:^ liminf ^ 12 p ,-^2 


(311.4) 
we get 


(311. b) 


and it remains to be shown that 


limsupF^:^ • (311.6) 

Case 1: Assume that liminfF'= iimsupF'= limF^= d < » ; integration 
of this relation F^^ d to F^dt yields by comparison with 
(311.4) d= 12u^'^ 9 and the proof is complete. 

Case 2: Assume that limF^= +c» and let d> 0 be arbitrarily large; 

then \iltimately F^> d and by integration F> dt , contra- 
dieting (311.4). 
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vlim inf F^< + 3 ^ 

Case 3; Assume finally that liminfF^< limsupF^, implying that 

There exist two positive numbers and dg and an 

infinite sequence of intervals t < t < t , both t and t eoine 

a B a B 

to zero, such that 


0 < 12^*2 < d, = FMt )< F'(t) < F'(t) = dg , t€ ]t ,t[ . (311.7) 

O' 8 O' 8 


Write 


F" 


4(- 


3 5 T^t^WS 


F'2 + 3^S-pN3-J[l_ 


and use (310.5) and (311.4) ; then 


|F"| < f- + F'^j 

and on It , t [ , since F'^ d~ : 

o 8 


Integration over [t yields 


V i. 

B B 

0 < dg - d, = F'(tg) - F'(0 = J F''‘(t)dt )F'' 


(t) 1 dt 


a 


a 


-t 


I t 

3 

iog~p = clog— ^ — 
a 8 


(observe that t < t < 0 ) ; this implies 
O' 8 


“ (V <^2 - d, 

log— — > , i. e. , ultimately 

“ 8 


-t 


O' 


-t 


s T > 1 


(311.8) 


From (311.7) it follows that 

t. 


F(tJ - F(t ) = j F'(t)dt > d- (t - t ) 

8 J ^ 3 a 

t 

a 


8 
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implying 


F(t ) 

gi 

-t 

6 


F(t ) 


-t 


+ > 


a 


-t 


+ d. 


a 


this last inequality is ultimately inconsistent, since 

both ^ , 


F(t ) F(t ) 

by (311.4) and “ 


-t 


B 


-t 


a 


by (311.7) 12g*® ^ , and 


-t 


by (311,8) ultimately — 1 . 

' B 

1 

Thus, J*® g* is true, and it follows that 

J*" ~ (y g*] ® t “ ® 


312 h^ remains bounded as t -* 0 : the energy h^ 
cluster is by (104.5) and (104.6) 

N 


of the 




k=l 


and by differentiation we get 


N 


& 


2^’"4,kri.kir.k * “i 


k=l 

substitute from the eqviation of motion (104.3): then 


+ u 


k=l 


k=l 


(311.9) 

i-th 

(312. 1) 


the first and ias<- term on the right cancel each other, thus 
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h' = m r ' P 



k=l 

and by integration 


N. 


h, = const 
I 


k=l 'J 


dt 


(312.2) 


Now J* 


. y \ Z Z S 

^ ^ \ ^ t ■ ^ implies V* < c t " ® 


and 


i£x.k' 


and since P is boiinded on [t, 0[ , we get 


0 


0 

If- 




= ct 


together with (312. 2) implying the boundedness of h . 

X 


313 The Liagrange -Jacobi equation for J 


from the definition (104.7) we get 

N. 


N, 


k=l k=l 


(312.3) 


(312. 4) 


and substituting V from (104.5) and for cITl. from the equation 

X ^X« K 

of motion (104. 3), 


N. 


N. 


j''‘ = 4V, - 2 ^ C„ ^ u) + 2 , r , 


k 


k=l 


k=l 


on the right, the third term is bounded, and the second term equals 

+2U (by Euler's theorem about homogeneous functions), and using 

IL 

h = V + U , where h is bounded (cf. 312), we find 
XXX X 
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J'' = -2U + b . (313. 1) 

Jo Jv 

As in (303, 3) we conclude (observing that J 0 by the assumption 

Xt 

of a collision singularity): 


n' 


= + b and U as t 0 imply: 

Xi jL 

+00 

ultimately does not change sign (and is negative) 
is ultimately mono tonicalfy decreasing to zero. 


(313.2) 


314 Estimate of (cf. 304) : 

(304. 1) , together with 

C, ,<cjf , (314.1) 

' JL,k JL 

yields 

|K^ I <cJf or equivalently , (314.2) 

' i X '^JL a ^ 


and integrating from t to t = 0 , observing that by (313.2) 
J = J^(t) = maxJ (t) , T € [t, 0[ , we have 

A> X/ Jo 


X 

K = K. + J^tb 
'**£;o j? •* 


eind 


i 

s 


- Kf + Jf tb 
I 1,0 SL 


(314. 3) 


315 ^ + d) : 

4 £ 1,0 1' H, 

analogous to the process used in 305 leading to (305. 1), we arrive at 


4-jy + s 2J V, 
A i I i I 


(315.1) 


by (313. 1): 2V = 2h - 2U - J'' + 2h - b ; using (314.3), 

Jo Jo SL Mt SL 

and observing that h is bounded, we get 

SL 




(315.2) 
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As in 307 for J* , it can be shown that (313. 1): j'' = -2U + b 

I i 

implies jf j" s c> 0 , hence 5 cJ j'' ; substitute this into the 

I I III 

last term of (315.2): 

^ ^ ^ '■'I-'-*" • 

2 

Finally, with j'' < J*'^+c , J^'^^ct'"^ (cf. (311.9)): 

(“t)j'' ^ c(-t)^ , and we get 

X* 

T-’r ^ v^r 

316 In full analogy to the development in 306 - 311, the following results 
can be derived: 


II 

0 

or 

= 0 ; 
l.o 

(316. 1) 

limj 2 

1 1 

II 

> 0 exists ; 

(316.2) 

\ ~(t“J 


9 

(316. 3) 

j; ~ U2,|j 


9 

(316. 4) 

t" /— ^ 

■'j ~l 3 % 

i 2 2 

3 t'3 

• 

(316. 5) 


317 Direct proof of K. = 0 : 

“*i;o «♦ 

we had seen that (cf. (312.3)) 

L.I < c(-t) ^ or equivalently i, ~ ^ » (317.1) 

and hy integration from t to t = 0 



results. Therefore 


K = / * m - r . X 

^ I Ji,ku,ji,k «*i. 


= t^b 


(317. 3) 


k=l 


i. e. , K -* 0 as t -> 0 

^ vvc 


(317.4) 


(317.4) together with (314.3) iraplies 
K =jftb and Kf = J ^ h 

leading to a simple derivation of the inequalV^y (315.4). 

318 Corollary. The estimate of K (cf. 314) can be dr eloped in 
full analogy with that of (cf. 304) : indeed, 

P, = 0 if C, =0 , k=1--N^ , (318.1) 

and writing 

Z«Z-C +c =w +C 

^X,x. vii, V '-‘Xk V 


and 


Iw +r |-^=w‘'^(l + r R ) , 


where R is bounded on ft ,t=^T , t^ sufficiently close to t^^ , 
Xhv ^ o 

(318.1) shows that all terms in P , not containing a factor C, 

^£,k ^ -jJ, V 

or C . cancel out. Since 

V 

1 
8 


Jl, V z 


it follows that 




and 




for the further steps see (304. 4) ff, . 
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319 In 319 - 323 we will show that each cluster's configuration is 
close to a central configuration if t < 0 is close to the instant 
t = 0 of a collision singularity; more precisely: the N bodies of 
the jJ-th cluster are close to a central configuration in the N^-body 
problem. 

We consider throughout the following an arbitrary, but fixed cluster, 
say the i-th. 


Let ilf be one of the components of C. , ; the conditions that the 
£-th cluster form a central configuration are 


= -iU 


I 9 \!( 


for all ijr 


(319.1) 


We introduce new dependent variables 



= t ^ r|t 


(320. 1) 
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'k = ‘’’Ci.k 


(320. 2) 


J~ = t " ^ 


(320, 3) 


U" = t ® U , 


(320.4) 


The conditions (319* 1) for a central configuration transform into 


for all i!(~ , 

the Lagrange “Jacobi equation ^ 

JL Ju At 

ej~'^ + -|-tj~' + J~ = -2U~ + b 


(320. 5) 


(320. 6) 


and the equation of motion (104.3b) into 

jJ, k\ 3 Uck 9 4*iky 

differentiating this last equation we get 


m ( t? cr"''' + ^ t c~' 

A,k ^k 3 -ik 9 «Mk , 




= -t-^ + t® m . ^(1 w ■*“ 4~ k * 

dt 9C, £, k dt w-ji, k 3 i, k*/»£,k 

4UXV 


321 From (316.3) - (316.5) it follows that 


4-1 - i 

tJ~'‘= - ®j+t , 

3 Ji i 


and similarly 


t 0 


(321.2) 



4l 


these relations together with (320, 6) imply that 

-yj~ ~ -U~ . (321.3) 

Substitute u ~ Cu'’ “T” ^ ^ the energy equation (312. 1); 

^SL f K 

we get 

" T«‘' " tJ" = 

implying by (321.2) and (321.3) that 

tC~'~0 . (321.6) 

By (316.5) or (316.3) and (321,3) U~ is bounded, therefore 
~ S U~ 

lim infer > 0 ; hence also — is bounded, and from (320,7) it 

It 

follows that 

Cr " i® bounded as t -♦ 0 , (321.6) 

4U-K 

In (320, 8) the first and second term on the right are bounded since 
tcr'‘ -*0 by (321.5), and the third term is obviously bounded; therefore 

t er"' i® hounded as t -• 0 , (321.7) 

322 (Tauberian) Lemma: If tf'(t) has a finite limit '9 and f''"(t) 

remains bounded as t -♦ 0 , then limt f (t) = -9 
as t -• 0 , (cf. Wintner [ 2 ] pp. 279 - 280) 

323 Applying 322 to (321.5), (321,6), and (321.7), we conclude that 

t cr" ~ 0 ; (323. 1) 

using this, together with (321.5) in the equation of motion (320,7), the 
asymptotic relation 
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2 5U~ 

— m , r ^ — 

9 i,k*»k 


or 


1 3J" 


9 ”^£,k' 


IZ 

3'li~ 


(323.2) 


results. The left side is — — and, multiplying both sides by 

-J~ ^ 0 and using (321,3), the desired result (cf. (320. 5)) 


r~ ^ 


-iu' 


MI 


(323. 3) 


follows. 
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400 Simultaneous binary collisions 


We shall consider the situation that there are pairs of bodies 

such that the two bodies in each pair collide at t = 0 , while the centers 
of mass of the pairs and the remaining Lq bodies approach distinct 
positions; thus 

Lj + Ijg = 1-1 = number of clusters, 

+ Lg = N = number of bodies. 


401 Change to a new coordinate system such that the motion of th^ 
second body in each cluster is referred to its first body. Introduce 
the new coordinate 


utit 2 1 


(401. 1) 


using m C, , + m^ _C, •, =0 (cf. (103.3)) , we find 
^ 1-A, 1 £,Zii£,,2 


and 


2 h ~ 1 1 

1.1^ = 1 ""il, Z^il, 2 


(401.2) 



402 The transformed equations of motion (104,3) read 






"'a. 2 

(402. 1) 


P. 



^ 4 / . 


the energy equation (104.6) changes into 


= V + 

I z 


m - m 
m , -r m ^ 

Z,\ l,Z 


m , *n , 

- V — ; (402.2) 


observing that 


m . m 

, 2 

I m 


2,1 2 ■' 


(402. 3) 


the Lagrange •Jacobi equation reads 


/ V' ^ 

(>:) -<7 • 


(402. 4) 


403 From the Lagrange -Jacobi equation (402.4) we derive now the 
a s y mpto tic expr ess ion 




( 1 


b) 


as t -4 0 : 


(403. 1) 


Abbreviate temporaril/ 


p2 = = 

JL 


(313.2) implies that H ultimately decreases monotonically to zero, 
and (316.3) that 


"2 = 


(403. 2) 


From the energy expression (402.2) follows " 0 , and since 

Ju ^ A/ ••^X/ 


-' = i p p' 0 by 


K -• 0 , we find 

— i - 


means 


f the identity 
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E is not identically zero since E is not constant. 

Ji S* 

Multiply the Lagrange -Jacobi equation 





+ b 


by 




and integrate from t = 0 to t, t<0 ; we get 


-\2 
~ I 


t 

= + f^lWmr)dr 

0 


(403. 3) 


Taking into account that V 

tiu.f u£i-9 Jucreases^ doe^ uci^ } 

ihe vo'Cae 

the last integral can be written as bH. ; hence, 

t 


-vs - 


=8r^.;u .b./t . 


and s ibstituting into the last term from (403,2) , we get 


H- = sr^./d . btl) 






yielding (403. 1) by integration. 
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404 In order to get an asymptotic expression for ^ as t — 0 , 

we first introduce a new independent variable s by 

1 

s = (404. 1) 

into the equations of motion (402. 1); denoting differentiation with 
respect to s by a prime ’ and substituting for ^ from (403. 1), 

J6 

we get 

-2sf +2E„ = s^b?+s®Q, . (404.2) 

Introduce further new dependent variables u^ by 

? - s u ; (404. 3] 

Ui ^ 


it follows that ^ = su , and since ^ = s^ b by (403. 1), both u^ 

a i ^ SL JL 

and u remain bounded and u can be written as u^ = sb • 

^ ***• 

The application of (404. 3) transform \ .04. 2) into 


u " = bu + Q 
'^JL 


(404. 4) 


substitute u = sb and integrate u" = sb + s^Q twice; then 

° ...jg — jg 


u=sa +s^b , 

««jg ^ 


(404. 5) 


P = s^a. + s^b 

**■ 


(404. 6) 


and by comparison with (403. 1) 




(404. 7) 
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405 Without restriction of generality, we can assume that all (finitely 
many) components of the are different from zero; this can always 

be achieved by a constant rotation of the coordinates. We now trans- 
form the equations of motion (104,2) and (402. 1) into a standard form 

of a system of diffex'ential equations of the first order, using again 

•> 

s = t*’ as the independent variable and introducing new dependent 
vnriables w , o). , W , and 0, by 


= 

;o 


+ s 


w 




= sW' 


(405. 1) 




S..J 



( 1 + SW ) 
Jo > J 


cu, = s w* 3=1,2,3 .(405.2) 
«*.je 


3 N ® 

Since Z = Z + s Z + s b , we have \v -• 0 as s -> 0 , 
***£ »»£;o '*‘^;o •>M.Z ^ 

and (404. 6) implies that w -* 0 , 
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Write for the further discussion (104.2) as 




(405. 3) 


and (402. 1) as 


^"=-r =~ + h 


(405. 4) 


where B and b are bounded as s,t-»0 and can be expanded \ 

■***£ I 

into power series of s, W . , and w . , which converge for 

ZlJ Z>) ! 

1 

sufficiently small absolute values of those quantities. 

406 Tromsformation of (405. 3) by (405. 1) gives 


so' = -Cl„ + 2W, + 9s" 


sW' = 0, 
j 2 v*Z 


407 From the first equation of (405. 2) we find 


= s^( a® + 2sx^ + s® b^ 


/ ’ 


where 


^ %; 2^;2 ^ 1 ; 3^£;3 ’ 


(406. 1) 


(407. 1) 


hence 


-a -6 

" "9 


H 


1 - 3a ® V + s^b 
Z ^Z 


■) 


(407. 2) 


Introduction of s and w into (405.4) and substitution from (407,2) 

^Z 

lead to 


-s 


s^w" . + 4sw' . = -2w. . + 6a, v« 4 sb* . , (407,3) 

j 2 ;j Z;j Z;j Z '"Z Z;j 
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where b-jf has the same properties as b. in 405. With uj 
(cf. {405, 2) ) we finally get the following system: 


scu' . = -3o). , - 2w. . + 6a. y. + sb*. 
^ fj £ ij JL Z Z Z tj 


j = l,2, 3 (407.4) 


S W' = U) . 


Z 9 } Z yj 


408 The full system of the transformed equations of motion (406, 1) 
and (407. 4) reads 


s n 


A;j 



+ 2W. . + s B* . 
Z‘,3 Zij 


sW' . 

Z 13 



j = l,2,3 
Z=l - Lj+La 


(408. la) 


-3 


s cu 


Zl3 


-3uj. . - 2w. . + 6a. Y. + s b* . 
Z‘,3 Z',3 Z Z Z',3 


s w' = tu. . ; 

Z 13 Z $3 


j = l,2,3 
Z = 1 — Li 


(408. lb) 


the nonlinear parts sB*_. and sb* . are power series in all 

£ 9J £ yJ 

variables, converging if their absolute values are sufficiently small . 


It is apparent that the matrix of the coefficients of the linear 
part on the right of (408. 1) is in "block form"; consequently, the set 
of its elementary divisors is the union of the sets of the elementary 
divisors of the individual blocks ( a block is the 6x6 matrix for a 
particular Z in (408. la) or (408, lb) ), 

The eigenvalues of a block in (408. la) are +1, +1, +1, -2, -2, -2 , 
and all elementary divisors are linear; similarly, the eigenvalues of 



50 


a block in (408. lb) are +1, -1, -1, -2. -2, -4 , and again all 

elementary divisors are linear. 

By a well-knowntheorem about systems of differential equations of the 
type (408. 1) it follows that the solution with the property that all W , 
n , w , and a' converge to zero as s 0 can be expanded into 
power series of s , which converge for sufficiently small |s| ; the 
solution depends on 4L«^ + SLg parameters. 

409 The theorem used in 408 has been proved by J. Horn: Ueber die 
Reihenentwicklung der Integrate eines Systems von Differ entialgleichun- 
gen in der Umgebung gewisser singul'drer Stellen. J. Reine Angew. 
Math. 116 (1896) 265 - 306 ; 117 (1897) 104 - 128 ; 254 - 266. For 
details and a survey of this and related problems see H. Dulac: 

Points singuliers des Equations differ entielles. Gauthier -Villars, 

Paris 1934. A modern simple proof (of a special case) can be found 
in Siegel [ 1 ] pp. 156 - 165. 
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410 Alternate proof of the boundedness of h . The principal tool 
for the study of the multiple binary collisions is the Lagrange -Jacobi 
equation (cf, (402.4)), enabling us to derive the asymptotic expression 

(403. 1) for 5 ; its applicability, in turn, is based on the boimdedness 

Jo 

of the energy h as t approaches the instant of collision. We deduced 

Jo 

this boundedness from general results abou ‘"^llisions (cf. 312). In the 

following we shall give a direct proof of the boundedness of h , 

Ju 

independent of the previous general results; for the method see 
Arenstorf [3]. 

According to 318, Q in equation (402. 1) can be written as 
5 b, provided that P sufficiently small, !• e. , that t is 

Jo *** Jo 

sufficiently close to the instant of collision t = 0. From 


h' = Q 


(cf. 312) 


(410. 1) 


= I l' b 


we get 


h'' ® 5 ( b*, to indicate j)ositive_ bounded function) 

Jw Jw 


and substituting for f from (402. 2) , 




and becaus e of 5 . 0 : 


h'; . 1 . ihj 


(for sufficiently small ? ) 

Jo 
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Denote by H (tj ) the maximum of ]h | on • il^en by 

integration from t^ to t and then to t^ , tQ < t ^ t^ : 

|h^(t) -h^(to)| ^ + H^(t,)^^(t - to)5 ^1 + H^(t,)^^(ti - to) . 

hence 

!h^(t)| ^ !\(to)| + - to) on [to.t,] 

and 

H^(ti ) = max |h^(t)| <: jh^(to)|+^l + H^{t^)^®(t^ - to) 

\ 1 

1 

For nonnegative a, b, x the inequality x a + b(l + x)^ implies 
a + ^b^ +b(l + a+ -^b^)^^a + b + ba^+b^; applying this, 
we find 

(\(to)l + ih -to)(l + |h_g(to)l^ + (ti - to)) . 

and it is obvious that the right side remains bounded as t^ -» 0 . 

411 We finally formulate the following 

Theorem: Let all variables be real and let the solution of the N-body 
problem be holomorphic on [t^, 0[ , ^ ^ there be 

at t = 0 a collision singularity such that several binary 
collisions (and only these) take place at this instant. Then 
the position vectors, distances, and the products 

t^ •(velocity vector) can be expanded into power series of 
1 

t^ . Therefore, the solution can be analytically continued 
in the real beyond t = 0 into t > 0 • 
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included; regularization of the equations of motion, non-newtonian 
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last ones in view of their importance for the collision problem. 
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